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By the first-principles electronic structure calculations, we find that the ground state of the Fe- 
vacancies ordered TlFei.5Se2 is a quasi- two- dimensional collinear antiferromagnetic semiconductor 
with an energy gap of 94 meV, in agreement with experimental measurements. This antiferromag- 
netic order is driven by the Se-bridged antiferromagnetic superexchange interactions between Fe 
moments. Similarly, we find that crystals AFei.5Se2 (A=K, Rb, or Cs) are also antiferromagnetic 
semiconductors but with a zero-gap semiconducting state or semimetallic state nearly degenerated 
with the ground states. Thus rich physical properties and phase diagrams are expected. 

PACS numbers: 74.70.Xa, 74.20.Pq, 74.20.Mn 



The discovery of high transition temperature super- 
conductivity in LaFeAsO by partial substitution of O 
with F atoms [1] stimulates great interest on iron pnic- 
tides. Other kinds of iron-based compounds were also 
reported to show superconductivity after doping or un- 
der high pressures [2|-[4[. A ubiquitous feature is that the 
parent compounds of these superconductors are antifer- 
romagnetic (AFM) semimetals|5j with either a collinear 
@,0| or bi-collinear (8l4Io| AFM order below a structural 
transition temperature. 

Very recently the superconductivity was discovered at 
about 30K in the potassium intercalated FeSe-layer com- 
pound Ko.8Fe2Se2 [U(. Soon after, the superconduc- 
tivity was also found in the Cs-intercalated compound 
Cso.8(FeSeo.98)2 [H| and (Tl,K)-intercalated compound 
(Tl,K)Fe^Se 2 [l3j]. These compounds have the ThCr 2 Si 2 
type structure (Fig. 1(a)), isostructural with 122-type 
iron pnictides BaFe 2 As2 [2j. But they should be regarded 
as a new kind of iron-based superconductors since they 
are chalcogenides rather than pnictides. We have done 
the first-principles electronic structure calculations for 
these materials [14|. We find that their parent com- 
pounds TlFe 2 Se2 and AFe 2 Se2 (A=K or Cs) are also 
AFM semimetals and the ground states are in a bi- 
collinear AFM order. 

An important feature revealed by the latest transport 
measurement is that the superconductivity in Tl and K 
intercalated FeSe materials is proximity to an AFM in- 
sulating phase [13j, similar as in high-T c cuprates. In 
particular, TlFe :E Se2 (1.3 < x < 1.7) is found to be an 
AFM insulator, with an activated transport gap of ~ 57 
meV for x = 1.5 [l3|- This has revived the discussion on 
the correlation effect in Fe-based superconductors. 

To clarify this issue, we have performed the first- 
principles electronic structure calculations on TlFei.5Se2 
and AFei. 5 Se 2 (A=K, Rb, or Cs). We find that the 
ground state of crystal TlFei.5Se2 is indeed an AFM 
semiconductor with a collinear AFM order (Fig. 4(b)) 
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FIG. 1: (Color online) TlFeccSe2 with the ZrCuSiAs-type 
structure: (a) a tetragonal unit cell containing two formula 
units without any vacancy; (b) schematic top view of the 
Fe-Fe square layer with one quarter Fe-vacancies ordered 
in rhombus (x = 1.5), in which there are two inequivalent 
Fe atoms, namely 2-Fe-neighbored and 3- Fe- neighbored Fe 
atoms, (c) schematic top view of the Fe-Fe square layer with 
one quarter Fe-vacancies ordered in square (x = 1.5). The 
filled circles denote the Fe atoms while the empty circles de- 
note the Fe-vacancies. The square and rectangle enclosed by 
the solid lines denote the unit cells. 



and an energy gap of 94 meV, and similarly AFei.5Se2 
is also an AFM semiconductor but with a zero-gap semi- 
conducting state or AFM semimetallic state close to the 
ground state. This is the first time theoretically to show 
there exists an AFM insulating state in Fe-based super- 
conductor materials. 

In our calculations the plane wave basis method was 
used [15]. We adopted the generalized gradient approx- 
imation (GGA) with Perdew-Burke-Ernzerhof formula 
[l6| for the exchange-correlation potentials. The ultrasoft 
pseudopotentials |l7j were used to model the electron-ion 
interactions. After the full convergency test, the kinetic 
energy cut-off and the charge density cut-off of the plane 
wave basis were chosen to be 800 eV and 6400 eV, re- 
spectively. The Gaussian broadening technique was used 
and a mesh ofl8xl8x9 k-points were sampled for the 
Brillouin-zone integration. In the calculations, the lattice 
parameters with the internal atomic coordinates were op- 
timized by the energy minimization. For TlFei.5Se2, the 
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FIG. 2: (Color online) Electronic structure of TlFei.sSe2 with 
the Fe- vacancies ordered in rhombus (Fig. lb) in the nonmag- 
netic state: (a) the band structure, (b) the Brillouin zone, and 
(c) the Fermi surface. The Fermi energy sets to zero. 




FIG. 3: (Color online) Electronic structure of TlFei.sSe2 with 
the Fe-vacancies ordered in square (Fig. lc) in the nonmag- 
netic sate: (a) the band structure, (b) the Brillouin zone, and 
(c) the Fermi surface. The Fermi energy sets to zero. 



optimized tetragonal lattice parameters are found in ex- 
cellent agreement with the experimental ones [13j . 

The early experiments suggested that the Fe-vacancies 
in TlFei. 5 Se2 are ordered in a rhombus structure as 
shown in Fig. 1(b) (lj). In order to test theoretically 
whether this structure is truly the ground state, we have 
also calculated the electronic structure of another pos- 
sible one-quarter Fe-vacancies ordered structure, namely 
the square-ordered vacancy structure as shown in Fig. 
1(c). In the nonmagnetic state, we find that the energy 
of the square-ordered vacancy structure is lower by about 
12 meV/Fe than the rhombus-ordered vacancy structure. 
It is noted that an Fe-vacancy only induces small local 
structural deviation (less than 0.04 A) from the origi- 
nal tetragonal one except for Tl atoms, which means the 
covalent bonding between Fe and Se atoms is rather ro- 
bust. Figs. 2 and 3 show the nonmagnetic band struc- 
tures and Fermi surfaces of TlFei.5Se2 in the Fe-vacancies 
rhombus- and square-ordered structures, respectively. In 
both cases, there are a number of electron and hole bands 
crossing the Fermi level. The Fermi surface nesting is not 
as strong as in the iron pnictides (H . 

However, in the true ground state which is in an AFM 
order, we find that the energy of the rhombus-ordered va- 
cancy structure is lower by 15.1 meV/Fe than the square- 
ordered one. Figs. 4 and 5 show a number of possi- 
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FIG. 4: (Color online) Schematic top view of four possible 
magnetic orders in the Fe-Fe square layer with one quarter 
Fe-vacancies ordered in rhombus: (a) checkboard Neel order 
in which the nearest neighboring Fe moments are anti-parallel 
ordered; (b) A-collinear AFM order in which the Fe moments 
are antiferromangtic ordered along the line without vacancies; 
(c) P-collinear AFM order in Fe moments are antiferromag- 
netic ordered along the lines with vacancies; (d) bi-collinear 
AFM order. The squares or rectangles enclosed by the solid 
lines denote the magnetic unit cells. 



ble magnetic orders in the Fe-vacancies rhombus- and 
square-ordered structures, respectively. Among these 
magnetic ordered states, we find that the Fe-vacancies 
rhombus-ordered structure with an A-collinear AFM or- 
der shown in Fig. 4(b) has the lowest energy. This 
rhombus order of vacancies agrees with the neutron 
measurement [18]. Our result further suggests that the 
rhombus-ordered vacancy structure is stabilized by an 
A-collinear AFM order. In the A-collinear AFM state 
(Fig. 4(b)), the Fe moments are antiferromagnetically 
ordered along the lines without Fe vacancies and ferro- 
magnetically ordered along the lines perpendicular. 

Similar as in the iron-pnictides [l9[ [5o[, we find that 
there is a small structural distortion in TlFei.5Se2. The 
lattice constant slightly expands along the AFM direc- 
tion and contracts along the ferromagnetic direction in 
the A-collinear AFM state. This leads to a small energy 
gain of ~1 meV/Fe. We also find that the Fe magnetic 
moments between the neighbor layers FeSe are antifer- 
romagnetically ordered. The energy difference between 
the ferromagnetic and AFM interlayer magnetic states is 
~4.9 meV/Fe. Thus the ground state of TlFei.5Se2 is an 
A-collinear AFM state with AFM interlayers of FeSe. 

Fig. 6(a) shows the ground state electronic band struc- 
ture of TlFei.5Se2- Unlike other parent compounds of 
iron pnictides or chalcogenides, TlFei.5Se2 is an AFM 
semiconductor, in agreement with experimental measure- 



FIG. 5: (Color online) Schematic top view of two possible 
magnetic orders in the Fe-Fe square layer: (a) c-AFM order in 
which the next nearest Fe moments are anti-parallel ordered; 
(b) h-AFM order in which half pairs of the next nearest Fe 
moments in anti-parallel while the other half pairs in parallel. 
The squares enclosed by the solid lines denote the magnetic 
unit cells. 
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FIG. 6: (Color online) Electronic band structure of TlFei.5Se2 
or AFei.5Se2 (A=K, Rb, or Cs) with one quarter Fe- vacancies 
ordered in rhombus (see Fig. 1(b)) and with an A-collinear 
AFM order (see Fig. 4 (b)) in the ground state, which 
is an antiferromagnetic semiconductor, (a) TlFei.sSe2; (b) 
KFei. 5 Se 2 ; (c) RbFei. 5 Se 2 ; (d) CsFei. 5 Se 2 . The Brillouin 
zone is shown in Fig. 2(b). Here the top of the valence band 
sets to zero. Note that T-N(r-N / ) corresponds to the anti- 
parallel(parallel)-aligned moment line. 



FIG. 7: (Color online) Projected density of states at the five 
Fe-3d orbitals around (a) 2-Fe-neighbored Fe atom and (b) 3- 
Fe-neighbored atom in the A-collinear AFM ordered ground 
state of TlFei.5Se2. Here the top of the valence band sets 
to zero. And x-axis is along the antiferromagnetic direction 
while 2/-axis is along the ferromagnetic direction in Fig. 4(b). 



merits. The energy band gap is found to be ~ 94 meV, 
which is also consistent with the gap value obtained by 
the transport measurement, 57.7 meV [l3|- The com- 
pound TlFei.5Se2 in the other magnetic orders is found 
mostly to be metallic. 

There are two inequivalent Fe atoms according to 
the number of neighboring Fe atoms if the Fe-vacancies 
are rhombus-ordered, namely 2-Fe-neighbored and 3-Fe- 
neighbored Fe atoms respectively (see Fig. 1(b)). By 
projecting the density of states onto the five 3d orbitals 
of Fe (Fig. 7), we find that the five up-spin orbitals 
are almost completely filled on both kinds of Fe atoms. 
This suggests that the crystal field splitting induced by 
Se atoms is small, similar as in iron-pnictides. The mag- 
netic moment formed around each Fe atom is found to 
be about 2.8/zb. This large Fe magnetic moments appar- 
ently results from the Hund's rule coupling ]j| 20|- The 



down spin orbitals are partially filled by d yz , d xy , and 
d x 2_ y 2 orbitals for the 2-Fe-neighbored Fe atoms or by 
d yz and d xy orbitals for the 3-Fe-neighbored Fe atoms. 
Such anisotropy among the five down spin orbitals in 
TlFei.5Se2 results from the Fe-vacancies. 

Inspection of the real space charge distribution in 
TlFei.5Se2 shows that there is a strong covalence bond 
between the neighbor Fe and Se atoms. This gives rise 
to an effective AFM superexchange interactions bridged 
by the Se atoms between the next nearest neighbor Fe 
moments, similar as in iron pnictides [HI, [lo|. If the 
energy of the nonmagnetic state is set to zero, we find 
that the energies of the ferromagnetic, checkboard Neel 
AFM, A-collinear AFM, and bi-collinear AFM states are 
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(-0.114, -0.217, -0.326, -0.270) eV/Fe for TlFei. 5 Se 2 . The 
magnetic moment around each Fe atom is found to be 
about 2.5-2.8/i#, varying weakly in the above four mag- 
netically ordered states, similar as in the iron-pnictides 
[l9l . l2o| . If we attribute the energy differences of these 
states entirely to the contribution of magnetic interac- 
tions between the spins S and model them by the sim- 
ple Heisenberg model with the nearest (Ji), next-nearest 
{J 2)1 and next-next nearest (J3) neighbor interactions 
and ignore the anisotropy caused by the Fe vacancies 
pjjj, we find that J x = 39 meV/S 2 , J 2 = 60 meV/S 2 , and 
J3 = 7 meV/S 2 for TlFei.5Se2- In obtaining these values, 
the contribution from itinerant electrons is ignored. Sim- 
ilarly we find that J x = 37 meV/S 2 , J 2 = 64 meV/S 2 , 
and J3 = 8 meV/S 2 for KFei.5Se2 (see below for more). 

To quantify the electronic correlation effect in 
TlFei.5Se2, we have performed a GGA+/7 calculation. 
We find that the energy band gap increases dramati- 
cally with the Coulomb interaction U . When U is larger 
than 2 eV, the energy band gap is over 190 meV, which 
is si gnifi cantly larger than the measurement value of 57 
meV[13]. This suggests that the magnetic ordering and 
the energy band gap are driven mainly by the exchange 
effect, rather than the correlation effect. 

Here we emphasize that the calculation convergency 
test needs to be elaborated. An enough high energy cut- 
off and a set of sufficient many k-points are required to 
ensure the correct ground state, namely an AFM semi- 
conductor obtained. Otherwise, the calculations always 
yield a metallic state rather than semiconducting state. 

We have further performed the first-principles elec- 
tronic structure calculations for AFei.sSe2 (A=K, Rb, 
or Cs). Like TlFei.5Se2, we find that crystal AFei.5Se2 
is an AFM semiconductor and its ground state is also A- 
collinear AFM state with AFM interlayers of FeSe. Figs. 
6 (b), (c), and (d) show the ground state electronic band 
structures for KFei.5Se2, RbFei.5Se2, and CsFei.5Se2, re- 
spectively. Accordingly, the energy band gaps are found 
to be 121, 69, and 26 meV for KFei. 5 Se 2 , RbFei. 5 Se 2 , 
and CsFei.5Se2, respectively. 

The calculations show that the energies of KFei.5Se 2 , 
RbFei. 5 Se2, and CsFei.5Se2 in the A-collinear AFM or- 
der but with the ferromagnetic interlayers of FeSe are 
higher than the AFM interlayers of FeSe by 8.9, 5.2, and 
2.8 meV/Fe, respectively. In such a case with the fer- 
romagnetic interlayers of FeSe, crystals TlFei.5Se2 and 
KFei.5Se2 both are still AFM semiconductors with indi- 
rect band gaps of 20 and 32 meV and direct gaps of 46 and 
62 meV, respectively. However, crystal RbFei.sSe2 be- 
comes a zero-gap semiconductor while crystal CsFei.sSe2 
becomes an AFM semimetal with the electron and hole 
carrier densities of 7.47 x 10 19 /cm 3 and 7.35 x 10 19 /cm 3 , 
respectively. These two states, zero-gap semiconducting 
and AFM semimetallic, will significantly influence the 
physical properties of RbFei.5Se 2 and CsFei. 5 Se2 since 
they are almost degenerated with the ground states. 



Moreover, we also find that the square-ordered vacancy 
structure with a c-AFM order (Fig. 5(a)) and AFM inter- 
layers of FeSe are metastable for crystal AFei.5Se2 (A=K, 
Rb, or Cs). Energetically, it is just higher than the 
ground state by 10.9, 7.9, and 5.0 meV/Fe for KFei.5Se 2 , 
RbFei.5Se2, and CsFei.5Se2, respectively. Such a square- 
ordered vacancy structure pattern should thus have a cer- 
tain probability of formation in the compound AFei.5Se 2 . 

In conclusion, we have performed the first principles 
calculations for the electronic structure and magnetic or- 
der of TlFei.5Se2. We find that the crystal structure 
of TlFei.5Se2 with the rhombus-ordered Fe vacancies in 
a collinear antiferromagnetic order is energetically the 
most stable. The ground state of TlFei.sSe2 is a quasi- 
two-dimensional antiferromagnetic semiconductor with 
an energy gap of 94 meV. Our result is consistent with 
both neutron and transport measurements. We further 
predict that crystal AFei.sSe2 (A=K, Rb, or Cs) is an an- 
tiferromagnetic semiconductor like TlFei.5Se2, and there 
is a zero-gap semiconducting state or antiferromagnetic 
semimetallic state nearly degenerated with the ground 
state for RbFei.sSe2 or CsFei.5Se2. 
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